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Abstract 

Molecular dynamics simulations were carried out in order to ascertain which of the poten- 
tial multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistant with 
its measured conductance. Estimates of the conductance obtained through counting ions that 
cross the channel and by solving the Nernst-Planck equation yield consistent results, indicat- 
ing that the motion of ions inside the channel can be satisfactorily described as diffusive.The 
calculated conductance of octameric channels is markedly higher than the conductance mea- 
sured in single channel recordings, whereas the tetramer appeal's to be non-conducting. The 
conductance of the hexamer was estimated to be 115±34 pS and 74+20 pS, at 150 mV and 
75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. 
On this basis we propose that the antiamoebin channel consists of six monomers. Its pore is 
large enough to accommodate K + and CD with their first solvation shells intact. The free 
energy barrier encountered by K + is only 2.2 keal/mol whereas Cl - encounters a substantially 
higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. 
Ion crossing events are shown to be uncorrelated and follow Poisson statistics. 
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Introduction 

Ion channels mediate and regulate transport of charged species across cell membranes. Not 
only do they play an essential role in the physiology of a cell but are also frequent drug targets. 
Despite their importance in biology and medicine, we know less about them than about nearly 
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any other major class of proteins. Only in the last decade were high-resolution structures of a 
number of ion channels solved through X-ray crystallography (see e.g. (1, 2)). However, a con- 
siderable gap still remains in understanding the structure-function relationship, as information 
revealed by X-ray crystallography is static and often incomplete. Model building and molec- 
ular dynamics (MD) simulations combined with X-ray structures can help close this gap, in 
principle, by providing insight into the dynamics of functional states and processes associated 
with ion conductance. 

To determine the relevance of MD simulations, they should be compared with experimental 
electrophysiological data, which directly measure the main function of ion channels — ionic 
conductance. At present, such comparisons arc very difficult to perform for large, eukaryotic 
ion channels or their bacterial homologs, as they require long simulations, most likely extend- 
ing to multi-microsecond time scales. Instead, we can test the approach and improve the com- 
putational tools using simple model channels. These channels not only inform us about how 
complex channels function, but also are of considerable interest in their own rights. Some arc 
viral channels which arc promising drug targets (3). Others, formed by non-genomic proteins 
from higher organisms, arc themselves therapeutic, as they exhibit antimicrobial activity (4, 5). 

In this study, we focus on ion channels formed by antiamoebin 1 (AAM), a member of a 
non-ribosomally synthesized family of fungal peptides called peptaibols (4, 6). Peptaibols arc 
short peptides that usually contain 15-20 residues and arc rich in non-standard amino acids, 
such as a-aminoisobutyric acid (Aib), D-isovaline (Iva) and hydroxyproline (Hyp). Many 
paptaibols form ion-conducting helical bundles that span cell membranes. Antiamoebin, first 
isolated from the fungal family Emericellopsis (7) and later found in other groups of fungi, is 
among the shortest peptaibols capable of forming channels, since it contains only 16 residues. 
The sequence of AAM is Ac-Phe-Aib-Aib-Aib-Iva-Gly-Leu-Aib-Aib-Hyp-Gln-Iva-Hyp-Aib- 
Hyp-Aib-Pro-Phol, where Phol stands for phenylalinol. None of the residues is charged at 
channel-forming, physiological conditions and, except for Gin and Hyp, all side chains arc 
relatively small and hydrophobic. The structure of a single AAM molecule, solved using both 
X-ray crystallography (8, 9) and NMR (10), is largely helical. 

The ion conductance properties of AAM have been studied in mixed palmitoyloleoylphos- 
phatidylcholine (POPC) / dioleoylphosphatidylcholine (DOPE) membranes (11). AAM chan- 
nels exhibit only one conductance level, which indicates a single structure. This is in contrast 
to alamethicin, which is the best studied peptaibol, both experimentally (12) and theoreti- 
cally (13, 14). It exhibits multiple conductance levels (15), which are thought to arise from 
a dynamic equilibrium between bundles containing different number of monomers (16). The 
microscopic conductance of AAM obtained from single channel recordings, adjusted to a 1 M 
concentration of KC1, is equal to 90 pS (11). It is believed that the channel exhibits some 
selectivity for K + over Cl“ (17). In addition to forming channels, AAM might also function 
as a carrier (9). 

The number of polypeptide chains needed to form conducting AAM channels is not well 
established. Measurements of macroscopic conductance as a function of peptide concentration 
led to the suggestion that the channel might be tetrameric (11). Under somewhat different 
assumptions about the equilibrium constants that govern channel formation, hexamers and oc- 
tamers are also possible. This led to a proposal that the pore of the tetrameric channel is 
too small to conduct ions and, instead, the channel is octameric (17). To resolve this issue, 
we built tetrameric, hexameric and octameric AAM channels, examined their short time-scale 
stabilities in the membrane, calculated their conductances and compared them with the con- 
ductance measured electrophysiologically. This approach, which directly links the structure of 
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the channel with its function, provides the means of determining the number of monomers in 
the conducting helical bundle. 

So far, channel conductances have been calculated mostly using methods that do not in- 
voke full atomistic details of the system, but instead rely on solving the Nernst-Planck (NP) 
or Poisson-Nernst-Planck (PNP and 3d-PNP) equations (18, 19) (see (20) for a recent review,) 
or on Brownian dynamics (19, 21, 22). Since important information about channel dynamics 
and energetics that influences conductance might be lost in these approaches, it is important 
to determine whether conductance can be accurately determined from all-atom simulations. 
In MD, channel conductance can be estimated directly by observing the ionic currents that 
develop in the presence of an applied electic field (23-26). Alternatively, conductance can be 
estimated from the electrodiffusion models (NP, PNP, 3d-PNP), in which it is assumed that 
ions move through the channel diffusively in response to the applied electric field and the un- 
derlying potential of mean force. In this study, we calculate conductance using a combination 
of both methods. Each method has its advantages and disadvantages. The direct method in- 
volves fewer approximations whereas the approach based on the Nernst-Planck equation can 
be applied to a broader range of conditions. Agreement between both methods indicates that 
the approximations underlying the diffusive model arc satisfactory for simple channels. 


Methods 

Initial Channel Models 

Initial structures of AAM bundles containing 4, 6 and 8 monomers were constructed using 
the m-Zdock program (27). For subsequent MD simulation we selected two octameric, one 
hexameric and one tetrameric channel using the procedure described in the Supplementary 
Material. One octameric model is a funnel-shaped bundle of parallel helices, as proposed 
by O’Reilly and Wallace (17). The lumen of the channel is hydrophilic due to the presence 
of the backbone carbonyl groups and the Gin 1 1 residues, which form a hydrogen-bonded 
ring pointing towards the center of the pore. The other candidate channel represents a class 
of well-packed structures in which the Gin residues are oriented towards the backbones of 
neighboring helices, forming a network of interhelical hydrogen bonds stabilizing the bundle. 
The selected hexameric structure is similar to that of O’Reilly and Wallace, but is tighter 
and has a smaller pore because it is built of a smaller number of monomers. The selected 
tetrameric channel is a compact structure in which the side chains of Gin 11 arc directed 
towards the neighboring helices. The minimum radius of the pore is 1.6 A, which is sufficient 
for transporting dehydrated, but not hydrated, K + ions. Structures in which the Gin residues 
are directed tow ards the center have pores that are too small to conduct ions. 

MD simulations. Several systems containing bundles of AAM embedded in a POPC bilayer 
placed between two lamellae of water were simulated using the program NAMD (28). Periodic 
boundary conditions were applied in three spatial directions. Long-ranged electrostatic forces 
were calculated using the particle mesh Ewald (PME) method with a grid of 128 x 128 x 
128. The r-RESPA multiple time stepping algorithm was used, in which the base time step 
was 1 fs and the long-range forces were updated every 4 steps. Water molecules were treated 
as rigid bodies using the SETTLE algorithm. Unless otherwise noted, the calculations were 
carried out in an NPT ensemble. The pressure was maintained at 1 atm by the Langevin piston 
with a period of 100 fs and a decay time of 50 fs. The temperature was kept at 300 K using 
Langevin damping with a damping coefficient of 1 ps 1 . Each system was initially equilibrated 
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using a multistep procedure described in the Supplementary Material. Then, 30, 20 and 45 ns 
equilibrium MD trajectories were obtained for the tetramer, hexamer and octamer, respectively. 

Potential energy functions. The CHARMM27 force field was used to describe the POPC 
phospholipids (29, 30) and water molecules were represented by the TIP3P model. CHARMM22 
was used to describe the standard amino acids in AAM (31). Non-bonded van der Waals pa- 
rameters and partial charges for the methyl groups attached to the o-carbon in Aib and Iva 
were taken from alanine, and the ethyl group in Iva was based on the parameters for the 13- 
carbon in valine and the methyl group in alanine. The parameters for the hydroxyl group in 
Hyp were taken from threonine. Finally, the C-terminal phenylalinol residue was built from 
phenylalanine with a -CH 2 OH blocking group on the a-carbon. 

To test the potential parameters for non-standard amino acids, a single AAM molecule 
was studied in hexane and methanol. The results were compared to the X-ray structure (9) 
and the NMR structure obtained in methanol solution (10). These simulations, described in 
the Supplementary Material, indicate that the potential functions used to describe the peptide, 
including those for atypical amino acids, arc suitable for modeling the AAM channel. 

Free Energy calculations. To compute the free energy profile (potential of mean force) en- 
countered by K + or Cl as it traverses the pore, an order parameter was chosen as the 21 - 
component (in the direction perpendicular to the membrane surface) of a vector joining a se- 
lected ion and the center-of-mass of the AAM bundle. The free energy profiles for both K + 
and Cl - along the order parameter were computed using the Adaptive Biasing Force (ABF) 
method (32-34), as implemented in NAMD (35). The method has been shown to be very effi- 
cient in calculating free energy changes along different dynamical variables (35, 36). Here, the 
transmembrane region was subdivided into several strata or ’’windows”, 5 to 8 A wide. Then, 
trajectories at least 10 ns long were obtained with an ion constrained in each window, and the 
free energy profile across the pore was constructed by integrating the average force over all 
windows. In addition, a 50 ns trajectory was computed in the absence of constraints. From 
this trajectory, equilibrium density profiles, p e q(z), were obtained for each ion. On both sides 
of the channel, but not in the central region, the densities arc sufficient to calculate the free 
energy profile, A(z) = kp,T In p eq (z), with satisfactory accuracy. Here, k is the Boltzmann 
constant and T is temperature. The total free energy profiles were obtained by combining the 
two methods for calculating A(z). ABF covered the range of 2 between -15 and 15 A, whereas 
density profiles were considered reliable in the ranges of [-20,12] and [12,20] A. 

Conductance Calculations. The single-channel conductance is defined as the ratio of the 
observed current, /, through the channel to the applied voltage (37), In MD simulations, the 
conductance can be computed directly by applying an external electric field to the system and 
observing the instantaneous current that develops (23, 26). 


I = 


1 

A tL z 
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where q, and y it) are are the charge and position at time t of atom i, respectively, and L z is the 
length of the simulation cell in the direction perpendicular to the membrane. The total current 
is obtained by integrating Eq 1. The current can also be estimated by counting the number 
of ions that cross the channel during the simulation. As shown in the Supplemental Material, 
these estimates yield the same conductance to within statistical errors. 

Simulations in the presence of a constant, external electric field, E ex t, were carried out in 
an NVT ensemble, as described by Aksimentiev and Schulten (23). The volume of the system 


4 



was fixed because fluctuations in the box length lead to fluctuation in the applied transmem- 
brane potential (23). Since the aqueous phase is conducting, the applied field corresponds to 
application of an external, transmembrane voltage V = E ext L z , where L z is the length of 
the simulation cell in the z-dircction. For the hexamer and the octamer, a 150 ns and 15 ns 
trajectories were generated, respectively, at a potential of 150 mV. This potential is higher than 
75 mV applied in the most reliable single channel recording experiments (11) in order to ob- 
serve a sufficient current for a meaningful estimate of the conductance and for determining 
whether ion crossing events arc uncorrelated. Slow relaxation of the channel was observed to 
occur over the first 50 ns of the trajectory, and this has been excluded from the analysis. 

Channel conductance can be also estimated by assuming that ions diffuse through the chan- 
nel in a one-dimensional potential of mean force and the applied field along z. If ion-ion 
correlations and the dependence of the diffusion coefficient, D, on z arc ignored, which will 
be shown to be good approximations, the flow of ions can be calculated from the NP equa- 
tion (37, 38). In the steady-state approximation, net flow of ions of a given type through the 
channel, j, is: 


j = -D 


( d p(z) Pp) dA (z) p(z ) dV (z) \ 
\ dz k B T dz k B T dz ) 


( 2 ) 


where p (z) is the linear number density of ions at position z in the channel (in units of number 
of ions/length) and V is the energy of the ion due to the applied voltage. Integrating this 
equation between z a and z b with the exponential integrating factor yields: 


. = D P (z a ) exp [(A (z a ) + V (z a ))/k B T\ - p (z b ) exp [(A (z b ) + V (z b ))/k B T] 

J f£ ex P[( A (z) + V{z))/k B T]dz 

In macroscopic calculations z a and zj, are set to the edges of the channel. Then it is assumed 
that A(z a ) = A(zb), and p(z 0 ) = p(z b ) is expressed in terms of the bulk concentration and the 
area in the mouth of the channel. (37, 38) The equation, however, is valid for any choice of z a 
and Zb inside the channel. Calculating j from this equation requires knowledge of A(z) from 
separate simulations in the absence of voltage. Once A(z) is available j can be calculated at 
different applied voltages. For two voltages across the membrane, A IV' and AV' J , the ratio of 
the corresponding currents, j 11 and j" is 

f_ _ IZ ex P [(^ ( z ) + yv ( z ))/ k BT ] dz 1 - exp (AV u /k B T) 
jV fZ exp [(Av (z) + W 4 (z))/&bT] dz 1 — exp (A V^/ksT) 

if z a and zt, arc at the edges of the channel. 

Alternatively, Eq. 2 can be integrated through dividing both sides by p (z) and integrating. 
This yields 


D 


] = 


P&b) 

E^r\p (v) 


In 


+ (A b (z) - A a (z ) + V b (z) - V a (z))/k B T 


(5) 


This equation is not practical for macroscopic simulations because the nonequilibrium den- 
sity profile, p(z), needs to be known, but it is useful for testing the validity of the diffusion 
equation against MD. If z a and zj t arc located at the edges of the channel all required informa- 
tion is obtained from a single MD simulation with applied voltage. We pursued both methods 
in our analysis. 
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Results 


Octameric and tetrameric channel models 

Simulations of the octameric channels, initiated from the structure proposed by O’Reilly 
and Wallace (17), yield a stable bundle both in the absence and the presence of an applied 
electric field. However, the conductance of this model estimated from Eq. 1 is 910 pS at 1 M 
salt concentration. A similar estimate is obtained from the NP equation using the approach 
described below. Since these estimates are an order of magnitude higher than the measured 
conductance it appeal's that the model does not correctly represent the AAM channel. The 
same is the case for the second octameric model considered here. Even though its initial 
structure had a smaller pore, and therefore was expected to exhibit lower conductance, the 
pore opened up and reached the same size as the pore in the first model during the first 10 ns of 
unrestrained simulation after equilibration. This yields a conductance that is also inconsistent 
with experiment. 

The tetrameric channel was simulated for 30 ns in the absence of an electric field. To 
prevent the channel from closing and allow for its structural adjustments to the presence of the 
ion, a K + ion was restrained in the middle of the pore for the first 9 ns. Once the restraints 
were removed the channel remained stable and retained a continuous water strand across the 
pore, but the ion was expelled from the channel. No crossing events were observed during this 
simulation or during the simulation in the presence of a 0.15 V applied potential. 

To increase the likelihood of generating an observable current, a higher potential of 0.5 V 
was applied. Assuming that the channel is Ohmic, a conductance of 90 pS should yield, on 
average, 14 crossing events in 20 ns under these conditions. In our simulations of this length 
the channel remained stable, but no crossing events were observed. This null result suggests 
that the tetrameric structure is non-conducting or its conductance is markedly lower than 90 ps. 

Structural stability and flexibility of the hexameric model 

In the reminder of this section we present results for the hexameric channel, shown in 
Fig. 1, which was simulated with and without electric field. In both types of simulations the 
system remained stable. The density profiles of its major components — water, phospholipids, 
protein and ions — as a function of the z-coordinate perpendicular to the membrane did not 
show any signs of drift. These profiles obtained from the simulation without electric field are 
shown in Fig. 2. The density profiles for water and POPC overlap significantly due to pen- 
etration of water molecules into the phospholipid head group region. Si mi lar profiles were 
obtained in simulations of a pure POPC-water system (39, 40). From the density profiles of 
the K + and Cl - it follows that the ions exhibit some charge layering at the water-membrane 
interface, with K + penetrating the membrane more readily than Cl - . This is primarily because 
negatively charged carbonyl and phosphate oxygen atoms, located closer to the hydrophobic 
core of the membrane than the positively charged NH 3 groups, interact favorably with positive 
charges, but not with negative charges. Recent experiments on solid-supported phosphatidyl- 
choline membranes showed a selectivity for the cations in KC1 and NaCl, although this depends 
on the particular ion pair (41). A smilar effect was noted in simulations of dipalmitoylphos- 
phatidylcholine (DPPC) membranes with NaCl (42). 

The length of the AAM bundle is only about 23 A (see Fig. 2), which is substantially 
less than the width of the hydrophobic core of the unperturbed membrane (~30 A). There is, 
therefore, hydrophobic mismatch between the peptide and the membrane. As predicted by 
the mattress model (43), this causes thinning of the membrane near the channel, as shown in 
Fig. S4 in the Supplementary Material. 
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The channels, which arc built of short helical monomers, appeal's to be quite flexible. This 
view is supported by the calculations of the pore radius along the channel axis, carried out with 
the aid of the HOLE program (44), shown in Fig. 3. The size of the pore undergoes substantial 
fluctuations, mostly due to inward and outward motions of the individual helices, with standard 
error varying between 0.5 and 1 A. Its narrowest part is located in the middle, near the Gin 1 1 
residues. This is a result of the overall shape of the channel and the bulkiness of Gin compared 
to the Aib and Iva residues that line other regions of the pore. At an applied voltage of 150 mV, 
the average minimum radius is initially equal to nearly 3.5 A. After 30 ns the channel narrows 
somewhat, to about 2.5 A, and the location of the minimum shifts slightly. In the absence of 
the electric field the pore is quite si mi lar. The Gin residues do not form a ring of hydrogen 
bonds that protrudes into the pore, as proposed by O’Reilly and Wallace. Instead, as detailed 
in the Supplementary Material, they form hydrogen bonds with water inside the pore or, less 
frequently, with the neighboring peptide. 

As shown in Fig. 2, the pore is sufficiently wide for ions to retain most of their first hy- 
dration shells inside the channel. Even in its narrowest region, only one water molecule is 
lost from these shells. For K + , this water molecule is replaced by the oxygen atoms from 
the side chains of Gin and Hyp residues with an additional, smaller contributions from back- 
bone oxygen atoms and the blocking alcohol of the phenylalinol. In the interfacial region, a 
slight decrease in the hydration of the K + ion equal, on average, to 0.5 water molecule is com- 
pensated by the carbonyl and phosphate oxygen atoms from the lipids. As a result, the total 
solvation number for K + remains fairly constant along the channel. The hydration number 
for Cl - also drops slightly in the membrane region. However, in contrast to K + , this drop is 
not compensated by polar groups in the protein lining the pore or in the lipids because these 
groups are also negatively charged. As a result, the solvation number decreases slightly, by 
approximately 0.5. 

Potential of mean force for ion transfer across the channel 

The potential of mean force, A(z), experienced by K + and Cl in the channel is shown in 
Fig. 4. This quantity has long been recognized as important for understanding the mechanism 
of ion transport (45, 46) and is a necessary ingredient of diffusion-based methods for estimating 
conductance. The free energy barrier associated with crossing the membrane through the chan- 
nel, equal to 2.2T0.5 kcal/mol and 5.0 ±0.7 kcal/mol for K + and Cl - , respectively, is markedly 
reduced compared to the barrier experienced by the ions in unassisted permeation (47). The 
difference in the barrier heights implies that the channel is selective for K + over Cl - . No local 
free energy minima are present for either ion, indicating the absence of binding sites inside the 
channel. Analysis of the trajectories provided no evidence for correlations in the motions of 
ions of the same or opposite charge. 

Ion Diffusion and Transport Properties 

Diffusion inside wide channels is usually well described by Fick’s law. However, if the 
channel radius is too small (~ 1 to 2 A) to allow permeants to pass each other inside the 
channel, their diffusion becomes single-file and Fick’s law no longer applies (48), as observed 
in computer simulations of ion transport through gramicidin, (46) the KcsA channel (49) or 
small sized carbon nanotubes (CNTs) (50). The pore in the hexameric AAM channel, which 
has a minimum radius of 2.5 A, appeal's to be sufficiently wide for ions not to move in 
a single file. This assertion was tested by analyzing the time evolution of the probability 
distribution that an ion located at position zq at time t = 0 is at position 2 at time t. This 
analysis, described in the Supplementary Material, reveals that at time scales longer than 10 ps 
ion diffusion is Fickian. The diffusion coefficient for K + near the center of the channel is 
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estimated at 0.5 x 10 -5 cm 2 /s. 

For Fickian diffusion, the diffusion coefficient, D z (z), can also be estimated from the force 
autocorrelation function of an ion that is held stationary at position z in the channel (51). In- 
side the channel the dependence of D z (z) on z is weak, of the order of the statistical errors. Its 
estimated values for K + and Cl arc 0.46 x 1 0 5 cnr/s and 0.70 x 1 0 5 cm 2 /s, respectively. 
For K + this estimate is quite close to the value obtained from measuring ionic displacement 
from the center of the channel, as described above. For comparison, the calculated diffusion 
coefficients in the bulk water region arc ~ 2.2, and 1.9 x 10 -5 cnr/s for K + , and Cl - , re- 
spectively and the corresponding, measured values of D for the same ions arc 1.96 and 2.03 x 
10 -5 cm 2 /s (37). Thus, the diffusion coefficient is reduced inside the hexameric AAM channel 
relative to its value in bulk water by approximately a factor of 3-5. A decrease in D for ions 
or water molecules was also observed in other channels, such as KcsA (49), several biological 
model channels (52) and gramicidin A (53). It might be caused by the increase in local friction 
due to the rough inside surface of a channel. Consistently with this view, D increases inside 
channels with smooth surfaces, such as CNTs (50). 

Channel Conductance 

A total flux of 10 K + ions was observed during the last 100 ns trajectory with an applied 
voltage of 150 mV and at 1 M KC1 concentration, which yields a channel conductance of 
109 pS. All crossing events occurred in the direction of the field. No Cl - ions crossed the 
pore during that time. This is expected considering that the free energy barrier to channel 
crossing is substantially higher for Cl - than for K + . The ion trajectories, from which the 
crossing events arc identified, are shown in Fig. 5a. The integrated current (see Eq. 1) is 
shown in Fig. 5b. The conductance estimated from the slope of a single straight line fitted 
to the integrated current (see Eq. 1), shown in Fig. 5b, is equal to 126 pS. Since ion crossing 
events appeal - to be uncorrelated they should be well described as a Poisson process. Then 
the probability distribution of waiting times (times between consecutive events) should decay 
exponentially with a characteristic constant. A, which defines the number of events in the 
total time interval. As seen in Fig. 5b, the cumulative distribution, P c (t), that waiting time 
between events was shorter than t follows the Poisson statistics very closely, indicating that 
ion transport through AAM can be described as a single, homogeneous Poisson process. The 
value of A determined through best fit is 10.5 ions/100 ns. The corresponding statistical error is 
3.24. This yields the conductance of 1 15±34 pS. A good agreement between different methods 
for calculating conductance suggests consistency and convergence of the calculations. 

Ionic current can be also calculated from the NP equation. Typically, this involves an area 
term, which is often identified with the cross-sectional area at the mouths of the pore (38). 
Information about the cross-sectional area of the rest of the pore is already included in the free 
energy profiles. Assigning the area term cannot be done unambiguously and, as a result, the 
accuracy of the calculated current might suffer considerably. To overcome these difficulties 
we used an alternative approach, in which only quantities calculated from the MD trajectory 
are used as input. These quantities are density profiles from equilibrium or nonequilibrium 
simulations. If the limits of integration in Eq. 3 or Eq. 5 are set to the edges of the channel 
than one of those two simulations is sufficient. For other limits they need to be combined. 
We explored both approaches. Eq. 5 yields the current of approximately 15 ions/100 ns and 
shows very good stability with respect to changes in the limits of integration (see Fig. S7 in the 
Supplementary Material). In contrast, estimates of j based on Eq. 3 appeal - to be quite sensitive 
to boundary effects and the detailed shape of A(z). As a consequence, we were unsuccessful in 
obtaining currents that exhibited comparable stability with respect to the limits of integration. 
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Using Eq. 4, we estimated that the current is reduced by a factor of 3 when the applied 
voltage is reduced from 150 mV to 75 mV. On the basis of the calculated ratio, which, in 
contrast to the absolute current, is fairly insensitive to the details of A(z), we estimated that 
the conductance at 75 mV is 74±20 pS, in good agreement with the experimental conductance 
of 90 pS. 


Discussion 

Relation between structure and conductance 

The number of monomers forming the channel can be estimated from measurements of 
macroscopic conductance as a function of protein concentration. Using a simple chemical 
equilibrium model, in which n units aggregate to form a conducting channel (54), the value of 
n for AAM was determined to be between 3 and 4 (1 1). If the aggregating units arc identified 
as monomers, this means that the channel is tetrameric (17), as trimers arc too small to con- 
duct ions. Alternatively, the formation of channels might be a two-step process, in which the 
initial aggregation of monomers into dimers is very strongly favored. This step is followed by 
self-assembly of dimers into conducting structures. If this mechanism is correct then the AAM 
channel should be a hexamer or an octamer. On this basis, O'Reilly and Wallace proposed 
an octameric model (17) arguing that the pore in a tetramer would be too small to conduct. 
Other aggregation mechanisms, governed by equilibrium constants that yield a more complex 
relation between n and the number of monomers in the channel, arc also possible. The interpre- 
tation of n is further complicated by possible coexistence of conducting and non-conducting 
aggregates. Thus, determining the number of monomers that form the channel solely on the 
basis of measuring macroscopic conductance might be unreliable. 

Octameric models yield conductances that arc markedly higher than those measured exper- 
imentally. Tetrameric aggregates, in turn, arc either too short to form transmembrane bundles, 
or have pores that arc too small to conduct ions, or arc too loosely bound to be stable. This 
conclusion is borne out by direct simulation of tetrameric bundles, which arc not observed to 
conduct ions, even in the presence of a large, 500 mV transmembrane potential. This conclu- 
sion may not apply to all peptaibols. In alamethicin, tetrameric structures have been associated 
with the low conducting level. This interpretation is supported by the measured conductance 
of tethered alamethicin tetramers (55). However, some theoretical studies cast doubt on the 
role of tetramers in the conductance of alamethicin (13, 14). Furthermore, even though the se- 
quences of antiamoebin and alamethicin arc related, their conductance properties arc, in many 
respects, different (56). Therefore, indirect evidence for conductance of tetrameric structures 
of alamethicin may not be relevant to antiamoebin. 

In contrast to tetrameric or octameric models, the hexameric structure yields satisfactory 
agreement between the calculated and measured microscopic conductance. At the applied volt- 
age of 75 mV and 1 M concentration of KC1 the conductance calculated using the NP equation 
is equal to 74 pS, compared to 90 pS measured in single channel recordings. Furthermore, 
the hexameric structure is, as expected, selective for K + . This is largely due to the favorable 
interactions of K + with oxygen atoms of Gin and peptide backbone atoms and to the deeper 
penetration of K + than Cl” into the membrane. On these basis, we propose that the AAM 
channel is built of six helices. 

Comparison of different methods for calculating conductance 
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Direct simulation of the conductance involves fewer assumptions than stochastic models, 
such as the NP equation. For pores that arc both large enough and stable enough to develop 
statistically significant currents at voltages of interest, this is the method of choice. Simulations 
on o-licmolysin, which has a much larger pore than AAM, were carried out at voltages between 
±120 mV to ± 1200m V (23). In this system good results for the current-voltage dependence 
were obtained in 10 ns simulations at each voltage. Simulations on OmpF (24) and OmpC (25) 
at an applied voltage of 1000 mV yielded good agreement with experimental conductances 
over a range of temperatures at lower salt concentrations (<5 M), but at higher concentrations 
the effects of temperature on conductance were underestimated. 

For low-conducting channels, such as AAM, obtaining statistically significant results from 
direct simulations in the most interesting range of applied voltages would require MD tra- 
jectories of several microseconds. In shorter simulations, the applied voltage must be larger 
than those usually used in experiments, sometimes even non-physiological. Mapping out the 
current-voltage curve, which is often available from electrophysiological measurements, for 
non-Ohmic, low-conducting channels using direct simulations becomes impractical. This lim- 
itation prevents stringent testing of theory against experiment. 

In contrast to the direct approach, the NP equation can be used to calculate the conductance 
at any desired voltage. It relies on an assumption that motion of the ion in the channel is 
diffusive. For channels with sufficiently large pores, such as AAM, this appeal's to be the case. 
It is also assumed that the free energy profile of the ions at potentials of interest can be obtained 
simply by adding the free energy of a linear voltage ramp to the equilibrium free energy profile. 
This means that the applied potential has only a small effect on the structure of the protein and 
membrane. The comparison of the pore radii (see Fig. 3) of the AAM simulated at 150 mV 
and without applied voltage indicates that this assumption appeal's to be valid. 

The limitations of continuum models have been discussed in the literature (37, 57, 58). In 
particular, the NP equation neglects ion correlations and the PNP or 3d-PNP equations are to 
be preferred (37, 58), although 3d-PNP calculations tend to overestimate conductance, due to 
a lack of dynamic ion-ion correlations (19, 46). In the present work, the computed free energy 
profile contains all interactions that an ion encounters in the channel in the absence of an exter- 
nal electric field, including ion-ion interactions in a mean-field approximation. As ions seem 
to cross the channel as uncorrelated events that follow Poisson statistics, it does not appear 
that contributions not captured in this approximation play a significant role in conduction and 
a more involved PNP approach is necessary. 

Even though the assumptions of the diffusion model are fulfilled, the agreement between 
the current calculated directly and estimated from the model is not as good as it might be 
expected. The discrepancies might be due to uncertainties in p(z). However, the consistency 
in j calculated at different integration limits argues against it. A more likely explanation is 
that the applied electric field, which should disappear only in the conducting phase, might 
extend beyond the range of the relatively short AAM channel. Then the voltage drop along the 
channel would be smaller than assumed and would lead to overestimating currents calculated 
from the NP equation. If the field extends from -25 to 25 A, which are close to the half-widths 
of most density profiles shown in Fig. 1, the current will be 12 ions/100 ns. The range of 
electric field can be in principle obtained from MD, but this would require markedly longer 
simulation times. 

The present simulations have been carried out using standard, non-polarizable models, 
although polarization effects have long been of concern (for recent reviews see (59) and (60)). 
For example, simulations of gramicidin with a polarizable force field yields a lower free energy 
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barrier to ion permeation and estimated conductances that are in qualitative agreement with 
experimental result (61). Polarization effects become increasingly important with the increase 
in the number of water molecules that arc removed from the first hydration shell and balanced 
by interactions with polarized protein atoms (62). In AAM, the ions permeate the channel with 
their first solvation shells essentially intact, thus effects of polarization on the computed free 
energy profiles are expected to be small. 


Conclusions 

We have investigated tetrameric, hexameric and octameric models of the antiamoebin trans- 
membrane channel in a water-POPC membrane system. The conductance of the hexameric 
model of AAM is estimated at 120-150 pS at 75 mV and 1 M KC1, which is in satisfactory 
agreement with the experimental value of 90 pS. The conductance observed in the octameric 
bundle is over 10 times larger, and the tetrameric model does not conduct even at voltages of 
500 mV. Consequently, we conclude that the conducting state of AAM is most likely hexam- 
eric. 

Ions traverse the pore with their solvation shells mostly intact, which in turn means that 
the free energy barriers to ions crossing the channel are small. The pore is sufficiently large 
to allow for Fickian rather than single-tile diffusion of ions. Other assumptions underlying 
the Nernst-Planck equation also appeal - to be fulfilled for simple channels. As a result, the 
estimates of ionic conductance obtained from the direct method that relies on counting ion 
crossing events and from the NP equation are in good agreement once the range of applied 
voltage is properly accounted for. This, in turn, allows for combining MD simulations at 
a single voltage and stochastic calculations to estimate conductance at any desired voltage 
without making explicit assumptions about the pore dimensions. 

Using computer simulations to link structural information about ion channels with electro- 
physiological measurements, which provide the most relevant functional data, appeal's to be 
a useful strategy for investigating simple ion channels. It provides a sensitive test of whether 
our understanding how channels are built and how they conduct ions is satisfactory. As com- 
putational capabilities expand, this strategy might become very valuable for studying more 
complex channels. 
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Figure Legends 

Figure 1: 

Side and top views of a model hexameric bundle of AAM. Each helix is colored separately, 
and the Gin 1 1 residues are show in CPK colors. 

Figure 2: 

Density profiles of water (dotted cyan), lipids (cyan), AAM (blue), K+ (red) and Cl— (green) 
ions. Shown above the density profiles are the average hydration and solvation numbers for 
a K + (red and dash-dotted red) and Cl~ (solid and dotted green) as a function of its position 
along the AAM6 channel. The hydration number is the number of oxygen atoms from water 
molecules that are located closer to the ion than the minimum between the first and second 
peaks in the ion-water oxygen radial distribution function. These distances are 3.4 A and 
3.75 A for K + and Cl - , respectively. Solvation number is defined in a similar fashion, but also 
includes oxygen and nitrogen atoms from the protein and oxygen atoms from the lipids that 
are within the same distances. 

Figure 3: 

Average radius of the channel, along the axis of the pore, for the AAM hexamer with (solid) and 
without (dotted) applied held. The average pore radius for the AAM octamer model without 
applied held (dashed) is shown for comparison. 

Figure 4: 

Free energy for transporting K + (solid) and Cl - (dashed) ions through the transmembrane 
AAM6 pore. 

Figure 5: 

(a) Trajectories of ions in the pore region for a 150 ns simulation with an external potential of 
0. 15 V. In the first 30 ns there were 4 Cl - trajectories (dotted) and 9 K + (solid) trajectories that 
cross the pore. During the remaining 120 ns simulation, after the channel relaxed to a structure 
with a smaller pore, 10 K + and no Cl - ions crossed the channel; (b) integrated ionic current as 
a function of time; (inset) The number of waiting times shorter than t predicted from Poisson 
statistics for the total time interval of 100 ns (solid line) and encountered in MD simulations 
(solid diamonds). 
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